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1,  Introduction 


The  fluorapatites  have  received  considerable  attention  recently  as  hosts  for 
rare-earth  ions,  and  transition-metal  ions  in  various  ionization  states 
[1-5].  In  the  fluorapatites  M5(B04)3F,  the  M  ion  can  be  Ca^'*'  or  Sr^"'',  which 
have  ionic  radii  matching  Ln^'^;  lanthanide  ions  can  replace  the  ion  [6]. 
(For  simplicity,  the  various  fluorapatites  M5(B04)3F  can  be  referred  to  as 
MBAP.)  Flowever,  the  replacement  of  the  ion  by  Lrfi'^  must  be  accom¬ 
panied  by  charge  compensation,  achieved  either  by  the  inclusion  of  some 
optically  neutral  impurity  ion,  such  as  Na"^  (2Ca^‘'’  Na"*"  +  or  Si'^"'' 

(Ca^"''  +  >  Si^"*"  +  or  by  vacancy  creation  (3Ca^'*'  +  a 

vacant  Ca^"*"  site).  Other  more  sophisticated  forms  of  charge  compensation 
can  be  and  have  been  suggested  [7-9].  In  the  analysis  reported  here,  I  con¬ 
centrate  on  two  types  of  charge  compensation,  which  were  proposed  by 
Maximova  and  coworkers  some  time  ago  [6,9].  I  also  consider  the  case 
where  the  charge  compensation  is  remote  and  unidentified,  and  the  case 
where  the  neodymium  ion  enters  the  M2  site  and  a  calcium  vacancy 
occurs  at  one  of  the  two  nearest  M2  sites  in  the  crystal. 

This  report  considers  the  fluorapatites  Ca5(P04)3F  (CPAP)  and  Sr5(P04)3F 
(SPAP),  referred  to  collectively  as  MPAP.  In  these  fluorapatites,  the  crystal 
belongs  to  the  P62,/m  (176)  space  group.  The  M  ions  occupy  two  sites:  Ml 
in  the  4/ site  with  C3  symmetry  and  M2  in  the  6h  site  with  Cg  symmetry 
[10].  The  latter  site,  M2,  is  the  subject  here,  and  I  assume  that  the  Nd^'*'  ions 
occupy  this  site,  which  is  the  site  that  is  generally  assumed  to  be  occupied 
by  the  rare-earth  ions  in  the  apatite  hosts. 

Three  of  the  six  M2  sites  (in  the  imit  cell)  of  the  P63/ m  space  group  are  (in 
hexagonal  crystal  coordinates)  (x,  y,  1/4),  (-y,  x-y,  1/4),  and  (y  -x,  -x, 
1/4),  and  the  three  others  are  given  by  (x,  y,  z)  — >  (-x,  -y,  -z);  the  crystal- 
field  components  are  always  computed  at  the  (x,  y,  1/4)  site.  The  group  Q 
has  only  one  operation,  which  is  a  reflection  in  a  plane  perpendicular  to 
the  principal  axis,  which,  in  this  particular  case,  is  the  c-axis  of  the  hexago¬ 
nal  crystal;  since  I  am  interested  only  in  the  symmetry-preserving  charge 
compensators,  I  restrict  my  attention  to  the  first  set  (z  =  1/4)  given  above. 
To  preserve  the  Cg  symmetry,  I  consider  the  charge  compensator  ions 
occupying  the  sites  with  z  =  1/4,  and  in  all  cases  I  take  the  origin,  the  loca¬ 
tion  of  the  Ln^'*'  ion,  at  (x,  y,  1/4).  Thus,  where  we  have  oxygen  compensa¬ 
tion,  I  assume  that  the  fluorine  ion  at  (0,  0, 1/4)  is  replaced  by  an  oxygen 
ion.  For  sodium  compensation,  I  assume  that  the  calcium  ions  at  either  (-y, 
X  -  y,  1  /4)  or  (y  -  X,  -X,  1  /4)  are  replaced  by  sodium  ions  (or  any  ion  with 
+1  values  that  might  occupy  the  site).  Both  these  t5^es  of  charge  compen¬ 
sation  are  discussed  in  some  detail  by  Maksimova  and  SoboT  [7].  The 
detailed  x-ray  data  on  MPAP  for  M  =  Ca  and  Sr  are  given  in  table  1  for  the 
representative  ions  in  each  site,  and  the  coordinates  (hexagonal)  for  the 
equivalent  sites  of  each  ion  are  given  in  the  International  Tables  [10]. 
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Table  1.  Ciystallographic  data  on  M^{BO^)^¥  (M  =  Ca,  Sr;  B  =  P,  V). 
Hexagonal  Cl}^  (PS^/nt)  176,  Z  =  2 


Ion 

Site 

Symmetry 

X 

y 

2 

t 

Ml 

4/ 

Q 

2/3 

1/3 

2 

2 

M2 

6h 

Q 

X 

y 

1/4 

2 

B 

6h 

Q 

X 

y 

1/4 

oi 

6h 

Q 

X 

y 

1/4 

Ro 

02 

eh 

Cs 

X 

y 

1/4 

Ro 

03 

Hi 

Q 

X 

y 

2 

Ro 

F 

2a 

CsH 

0 

0 

1/4 

-1 

X-ray  data  on  M5(B04)3F 

M 

B 

a  {A) 

c(A) 

Zmi 

^M2 

^M2 

Xb 

Xoi 

Toi 

Cab 

P 

9.3973 

6.872 

0.0010 

0.00712 

0.24827 

0.36895 

0.3985 

0.4849 

0.3272 

Sr= 

V 

10.0077 

7.4342 

0.0004 

0.24943 

0.01059 

0.3682 

0.3982 

0.4706 

0.5962 

Sr<: 

P 

9.7153 

7.2810 

-0.0002 

0.23910 

-0.01445 

0.3992 

0.3685 

0.3306 

0.4810 

M 

B 

Xo2 

^02 

Xo3 

^03 

^03 

Cab 

P 

0.4667 

0.5875 

0.2575 

0.3421 

0.0705 

Sr^ 

V 

0.4835 

0.3160 

0.2498 

0.3403 

0.0671 

Sr^^ 

P 

0.5356 

0.1183 

0.3441 

0.2612 

0.0788 

^Two  values  of  qo  have  been  chosen,  -2  and  -1.8,  with  qp  = -3-  Aqo  in  each  case. 

M.  Hughes,  M.  Cameron,  and  K.  D.  Crowley,  Structural  Variations  in  Natural  F,  OH,  Cl  Apatites,  Am. 
Mineral.  74  (1989),  870. 

Private  communication,  R.  E.  Peale,  University  of  Central  Florida. 

2.  Theory 


In  all  the  calculations  reported  here,  the  crystal-field  Hamiltonian  for  the 
4/^  configuration  is  given  by 

n  ^  N 

^  ^  ^nm  /  (1) 

n  even  m  =  -n  i  =  l 

where  the  are  the  crystal-field  parameters  and  the  C„^  are 

^nm(^  ~  (2n  +  1)  (Pi)  ,  (2) 

with 


and  the  are  the  ordinary  spherical  harmonics.  In  Cg  symmetry,  the  val¬ 
ues  of  m  are  restricted  by  n  -i-  m  =  0,  ±2,  ±4,  with  \m  \  <n  and  n  =  2, 4,  and 
b-  The  parameters  entering  into  the  free-ion  Hamiltonian  are  from  aqueous 
solution  [11],  and  their  values  as  well  as  the  details  of  the  computation  are 
given  elsewhere  [12,13].  The  crystal-field  splitting  is  predominantly  deter¬ 
mined  by  the  n-even  from  the  discussion  above,  the  number  of 
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for  Cj  symmetry  is  limited  to  14,  with  B22  real  and  positive;  these 
properties  of  B22  ^re  simply  achieved  by  a  rotation  about  the  z-axis  (c-axis). 
In  the  three-parameter  crystal-field  theory  [14],  the  crystal-field  param¬ 
eters,  Bnf„,  are  related  to  the  crystal-field  components,  by 


^nm  ~  Pn^n 


(3) 


where  the  p„  values  represent  the  effective  values  of  (r”)  for  the  rare-earth 
ion.  Values  of  p„  have  been  tabulated  for  the  entire  rare-earth  series  [12]. 
The  crystal-field  components  are  given  by 


where  is  the  effective  charge  in  units  of  electronic  charge  on  the  ion  at  Ry 
[15],  and  the  sum  is  over  all  the  ions  in  the  solid.  Equation  (4)  has  been 
used  to  calculate  the  for  the  Ca2  site  in  CPAP  and  the  Sr2  site  in  SPAP 
with  the  x-ray  data  reported  in  table  1.  Two  values  of  the  oxygen  charge 
have  been  used:  qQ  =  -2  (for  which  -  2  and  qp  =  5)  and  qQ  =  -1.8  (for 
which  qo,  =  2  and  qp  =  -4.2),  which  in  both  cases  maintain  the  total  charge 
of  the  PO4  molecule  at  -3  {qp  =  -3  -  4qo)-  The  value  of  -1.8  for  the  oxygen 
charge  was  chosen  because  this  value  best  fit  the  experimental  value  of  the 
point  charge  determined  experimentally  for  rare-earth  ions  in 
Y3AI5O12  (yttrium  aluminum  garnet — ^YAG).  Table  2  has  the  results  of  the 
calculations  of  the  for  the  Q  sites  at  {x,  y,  1/4)  in  CPAP  and  SPAP  with 
no  charge  compensation.  The  for  CPAP  given  in  table  2  should  corre¬ 
spond  to  the  conjecture  of  Maksimova  and  SoboT  [7]  for  low  concentration 
of  Nd^"*"  ions. 

To  calculate  the  effect  of  replacing  the  nearest  F“  site  at  (0, 0, 1/4)  by  an 
site,  I  calculate  the  A^m  for  a  negative  charge  at  this  site  and  add  these 
to  the  corresponding  values  given  in  table  2. 

Table  2.  Even-«  crystal-field  components,  (cm"VA"),  for  M2  site  at  (x,  y,  1/4)  in  CPAP  and  SPAP 


(no  charge  compensation).^ 


A 

^7tm 

CPAP” 

SPAP^ 

^0  - 

-1.8 

^0- 

-2 

4o  = 

-1.8 

^0  = 

-2 

Re 

Im 

Re 

Im 

Re 

Im 

Re 

Im 

^20 

9481 

— 

11398 

— 

6535 

— 

8081 

— 

^22 

-2048 

-506.5 

-2384 

-717.6 

527.2 

-1664 

483.1 

-2034 

^40 

2031 

— 

2012 

— 

1902 

-1911 

— 

— 

^42 

2829 

-1801 

3187 

-1983 

-1815 

1053 

-2021 

1204 

A44 

589.8 

-2164 

369.5 

-2417 

1124 

-1298 

1386 

-1251 

^60 

90.85 

— 

121.2 

— 

113.2 

— 

141.5 

— 

^62 

-329.9 

-284.8 

-355.9 

-315.5 

-3.465 

-321.1 

-7.119 

-349.6 

^64 

156.0 

70.42 

193.1 

79.35 

-113.7 

-39.64 

-135.9 

-55.88 

^66 

64.57 

-158.4 

85.17 

-175.0 

47.50 

179.5 

63.08 

197.7 

and  y  for  Ca2  and  Sr2  sites  are  given  in  table  1. 

’’Smallest  Ca2  -  F  distance  =  2.3109  A,  smallest  Ca2  -  Ca2  distance  =  4.0025  A  (two  ions). 
'’Smallest  Sr2  -  F  distance  =  2.3962  A,  smallest  Sr2  -  Sr2  distance  =  4.1504  A  (two  ions). 


7 


Consider  for  a  moment  CPAP:  the  nearest  F~  to  the  Ca2  site  is  at  2.3109  A, 
and  the  next  nearest  F~  ions  are  at  4.1434  A;  these  facts  would  lead  one  to 
surmise  that  most  of  the  contribution  of  all  the  fluorine  ions  is  approxi¬ 
mately  due  to  this  nearest  ion.  That  this  is  the  case  is  shown  in  table  3  for 
CPAP;  if  a  rough  estimate  of  the  effect  of  the  nearest  F“  being  replaced  by 
an  ion  is  needed,  one  can  obtain  it  by  adding  in  the  contribution  of  all 
the  F“  ions  to  the  total.  That  is,  simply  double  the  contribution  of  the  F" 
ions  to  the  crystal-field  components  However,  I  do  not  make  that 
approximation  here  but  evaluate  the  crystal-field  components  for  a  single 
negative  charge  at  (0,  0, 1/4)  and  add  this  to  the  total  for  the  site  at  (x, 
y,i/4). 


Table  3.  Crystal-field 
components  for  a 
single  ion  at  (0, 0, 

^nm 

Single  F“ 

AUF- 

ions 

Re 

Im 

Re 

Im 

1/4)  and  total  F"  ion 

An 

-386 

15374 

^21 

17646 

contribution  to 

^20 

^706 

-2972 

crystal-field 

2^22 

-5756 

-289 

^454 

-318 

components  at  Ca2 

^31 

44.2 

-1763 

31.3 

-1247 

site  at  (x,  y,  1/4)  in 

^33 

171 

-2270 

179 

-2327 

Ca5(P04)3F  (CPAP).= 

^40 

661 

644 

2^42 

696 

34.9 

604 

30.3 

A44 

917 

92.3 

934 

92.9 

^51 

-6.55 

261 

-6.86 

274 

^53 

-21.2 

281 

-20.3 

266 

^55 

^7.3 

375 

^7.5 

375 

^60 

-103 

— 

-108 

— 

-106 

-5.30 

-no 

-5.54 

^64 

-115 

-11.6 

-112 

-11.4 

^66 

-155 

-23.5 

-155 

-23.5 

All 

1.05 

^1.7 

1.06 

-42.3 

^73 

3.26 

-43.3 

3.35 

-44.6 

^75 

6.00 

^7.6 

5.95 

-47.2 

An 

11.4 

-64.4 

11.4 

-64.4 

^AllA„„ 

values  in  uruts  of  cm' 

-VA”. 

3.  Calculation  of  Crystal-Field  Parameters,  Bnm,  for 
Various  Types  of  Charge  Compensation 

The  types  of  charge  compensation  considered  are  given  in  table  4,  and  the 
labels  given  in  table  4  that  identify  the  type  of  charge  compensation  are 
used  hereafter.  As  an  example  of  the  method  of  obtaining  crystal-field 
parameters  for  Nd3+,  I  give  the  details  of  the  calculation  for  the  Cl  type  of 
charge  compensation  in  table  4.  The  results  given  in  table  2  are  for  a  crystal 
where  each  ion  in  the  solid  is  in  the  position  given  in  table  1.  Then  to 
account  for  an  Na"*"  ion  replacing  a  Ca^'*'  ion  at  (— y,  x  —  y,  1/4),  we  can  calcu¬ 
late  the  Anm  for  a  single  charge  of-1  at  that  point  and  add  these  to  the 
appropriate  of  table  2.  After  this,  the  resulting  are  rotated  so  that 
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A22  is  real  and  positive  [16].  Using  the  rotated  we  can  calculate  the 
crystal-field  parameters  for  Nd^'*'  by  using  equation  (3).  For  the  D  types  of 
table  4,  we  can  create  a  vacancy  at  the  two  Ca^"''  sites  by  choosing  a  -2 
charge  at  the  D1  or  D2  site  in  the  table,  and  the  resulting  are  treated  as 
in  the  Cl  case  discussed  above.  The  results  for  all  the  cases  of  table  4  are 
given  in  tables  5  and  6.  The  invariants,  I„(B),  given  by  [17] 


m= 


y  B*  B 

m  =  -n 


1/2 


(5) 


were  calculated  and  are  also  given  in  tables  5  and  6.  The  results  given  in 
tables  5  and  6  show  that  the  crystal-field  parameters,  for  the  different 
types  of  charge  compensation  vary  considerably;  except  for  the  B  type  (F^ 
however,  these  variations  have  little  effect  on  the  I„(B)  for  the  dif¬ 
ferent  forms  of  charge  compensation.  The  result  found  here  for  six  of  the 
possibly  numerous  charge  compensation  possibilities  is  indicative  of  the 


Table  4.  Types  of 
charge  compensation 
considered  in  text.^ 


Type  Remark 

A  No  charge  compensation 

B  Ca^'^  and  F"  are  replaced  by  at  (0, 0, 1/4),  and  Nd^'*'  at  (x,  y,  1/4) 

Cl  Ca2+  at  (“-y,  x-y,\  /4)  replaced  by  Na-*- 

C2  Ca^'*'  at  (y  -  X,  -x,  1/4)  replaced  by  Na"^ 

D1  Ca^"^  at  (-y,  x  -  y,  1  /4)  removed  (a  charge  -2  at  this  site) 

D2  Ca^"^  at  (y  -  x,  -x,  1/4)  removed  (a  charge  -2  at  this  site) 


^Crystal  field  is  calculated  for  the  site  at  (x,  y,  1/4),  with  x  and  y  for  the  Ca2  site 
given  in  table  1.  In  aU  cases  read  Ca^'^  as  meaning  Ca^*^  or  Sr^*^.  (Na^^  can  be  taken  as 
any  possible  replacement  ion  with  valence  of  +1.) 


Table  5.  Nd^"^  crystal-field  parameters,  (cm“^),  and  rotation  invariants,  IJS)  (cm  ^),  for  different 
types  of  charge  compensation  in  CPAP  with  =  -1.8  and  -2.^ 


^nm 

^0  - 

-1.8 

II 

0 

A 

B 

Cl 

C2 

D1 

D2 

A 

B 

Cl 

C2 

D1 

D2 

^20 

1617 

815 

1463 

1463 

1309 

1309 

1944 

1142 

1790 

1790 

1636 

1636 

B22 

360 

1338 

458 

505 

672 

601 

425 

1399 

511 

572 

642 

739 

^40 

1173 

1555 

1197 

1197 

1222 

1222 

1162 

1544 

1186 

1186 

1211 

1211 

ReB42 

-1337 

-1922 

-1795 

-907 

-611 

-1951 

-1432 

-2087 

-1933 

-1025 

-2141 

-778 

ImB42 

1402 

1221 

788 

1709 

1816 

360 

1627 

1392 

1032 

1908 

564 

2020 

ReB44 

-281 

611 

711 

-908 

-1106 

1152 

-593 

394 

401 

-1104 

996 

-1266 

ImBqq 

-1264 

-1348 

-1115 

-879 

-541 

-701 

-1282 

-1483 

-1384 

-831 

-1081 

^75 

^60 

144 

-19.5 

141 

141 

137 

137 

193 

28.7 

189 

189 

186 

186 

ReB62 

618 

735 

448 

684 

699 

316 

686 

791 

531 

746 

393 

761 

ImB62 

314 

389 

527 

128 

24.3 

616 

317 

416 

537 

146 

645 

44,5 

ReB« 

272 

82.4 

201 

217 

156 

92.1 

326 

146 

286 

258 

175 

188 

ImB64 

-16.5 

78.5 

188 

-165 

-225 

262 

-64.3 

74.3 

172 

-211 

288 

-277 

ReBee 

90.7 

224 

-211 

253 

274 

-274 

127 

217 

-197 

281 

-308 

312 

ImB«, 

256 

233 

173 

104 

-13.2 

-10.9 

282 

254 

242 

133 

49.3 

14.1 

hm 

1696 

2061 

1600 

1628 

1617 

1560 

2035 

2285 

1930 

1964 

1871 

1941 

/4(B) 

3498 

4143 

3552 

3480 

3464 

3606 

3839 

4437 

3894 

3822 

3950 

3807 

hiB) 

1130 

1272 

1131 

1135 

1139 

1132 

1262 

1369 

1263 

1267 

1264 

1272 

^See  table  4  for  explanation  of  symbols  used  for  different  types  of  charge  compensation. 
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Table  6.  crystal-field  parameters  (cm“^),  and  rotation  invariants,  I„(B)  (cm"^)  for  different 
types  of  charge  compensation  in  SPAP  with  qQ  =  -1.8  and  -2.^ 


qo  -  -1.8 

‘)o  =  - 

2 

^nm 

A 

B 

Cl 

C2 

D1 

D2 

A 

B 

Cl 

C2 

D1 

D2 

^20 

1115 

395 

976 

976 

838 

838 

1379 

659 

1240 

1240 

1102 

1102 

B22 

298 

1178 

422 

397 

570 

533 

357 

1234 

485 

443 

633 

567 

^40 

1099 

1417 

1119 

1119 

1139 

1139 

1104 

1422 

1124 

1124 

1144 

1144 

R6B42 

-896 

-1299 

-602 

-1160 

^36 

-1241 

-946 

-1404 

^72 

-1249 

-508 

-1367 

ImB42 

-816 

-745 

-1049 

-413 

-1127 

-160 

-975 

-857 

-1179 

-582 

-1257 

-295 

R6B44 

-99.2 

515 

-672 

645 

-846 

932 

-390 

333 

-839 

387 

-971 

813 

ImB44 

987 

1015 

693 

783 

408 

457 

1005 

1120 

640 

1033 

341 

781 

^60 

180 

52.8 

177 

177 

175 

175 

225 

97.7 

222 

222 

219 

220 

ReB62 

485 

600 

513 

386 

507 

305 

538 

644 

558 

453 

554 

370 

ImB62 

-159 

-197 

4.26 

-335 

89.7 

-410 

-139 

-204 

6.52 

-323 

88.5 

-416 

ReB64 

184 

29.8 

177 

89.3 

140 

5.55 

233 

77.8 

209 

160 

166 

65.8 

ImB64 

-52.6 

-91.7 

73.8 

-172 

131 

-197 

-17.8 

-94.6 

106 

-173 

165 

-230 

113 

232 

289 

-218 

283 

-294 

176 

243 

322 

-176 

322 

-315 

IrnBf^ 

-273 

-276 

-60,3 

-199 

84.2 

-18.2 

-279 

-294 

-71.0 

-279 

72,0 

-99.6 

m 

1192 

1712 

1144 

1126 

1162 

1127 

1469 

1865 

1417 

1389 

1419 

1363 

14(B) 

2472 

2355 

2458 

2519 

2446 

2565 

2690 

3188 

2678 

2738 

2665 

2786 

k(B) 

895 

735 

897 

896 

899 

897 

997 

1115 

1000 

999 

1002 

1000 

“See  table  4  for  explanation  of  symbols  used  for  different  types  of  charge  compensation. 


number  of  lines  that  appear  in  the  absorption  and  emission  spectra,  as  was 
for  example  observed  by  Zounani  et  al  [18]  for  in  SPAP. 

The  crystal-field  parameters  given  in  tables  5  and  6  were  used  in  calculat¬ 
ing  the  energy  levels  for  Nd^"*".  The  details  of  the  procedure  used  in  the  cal¬ 
culation  are  given  elsewhere  [19].  The  calculation  covered  the  lowest  13 
multiplets,  with  a  total  of  60  energy  levels.  Instead  of  presenting  all  the 
results,  which  are  quite  extensive,  I  have  chosen  to  give  only  the  width  in 
energy  of  each  multiplet  where  such  a  quantity  can  be  reasonably  deter¬ 
mined.  The  results  for  all  the  given  in  tables  5  and  6  are  presented  in 
tables  7  and  8.  The  widths  of  the  Hj  (J  =  9/1  through  15/2)  and  the  ^^3/2  are 
well  determined,  and  greater  significance  should  be  given  to  these  multi¬ 
plets  than  to  the  widths  of  the  higher  multiplets  when  one  compares  these 
values  to  experiment.  Because  of  the  very  large  magnitude  of  the  B„^,  in 
particular  B2ni,  the  higher  multiplets  are  severely  mixed.  In  a  number  of 
cases,  the  mixture  was  so  great  that  the  number  of  levels  of  a  multiplet 
allowed  by  group  theory  was  exceeded.  This  occurred  numerous  times  for 
the  ^83 /2  levels,  which  had  three  levels  identified  (the  levels  are  identified 
by  the  largest  coefficient  of  the  state  in  the  wave  function  for  that  state), 
where  by  group  theory  there  should  be  only  two.  This  result  does  not  indi¬ 
cate  an  error  in  the  calculation;  it  simply  means  that  the  crystal  field  is  so 
strong  that  in  these  spectral  regions  it  renders  meaningless  the  method 
used  to  identify  the  energy  levels  by  assignment  to  particular  multiplets  in 
the  weak  crystal-field  limit.  Since  the  crystal-field  parameters  given  in 
tables  5  and  6  are  reasonable  approximations  to  the  values  that  would  be 
obtained  by  fitting  the  spectra,  these  results  indicate  the  difficulties  to  be 
expected  in  the  analysis  of  experimental  data. 
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Table  7.  Total  crystal-field  splittings,  AEj  (cm”^),  of  each  /  multiplet  of  Nd^'*'  in  CPAP  in  site 
with  qQ  =  -1.8  and  -2.® 


^0 

=  -1.8 

qo  = 

-2 

[(L,S)/] 

A 

B 

Cl 

C2 

D1 

D2 

A 

B 

Cl 

C2 

D1 

D2 

%/2 

649 

747 

632 

645 

645 

624 

732 

795 

714 

729 

702 

729 

%l/2 

506 

608 

482 

504 

503 

460 

581 

632 

557 

578 

535 

577 

%3/2 

589 

714 

564 

586 

584 

542 

676 

746 

652 

673 

630 

671 

%5/2 

849 

994 

823 

850 

849 

800 

980 

1071 

955 

980 

931 

981 

‘*^3/2 

228 

394 

235 

238 

258 

254 

262 

439 

270 

276 

289 

298 

^53/2*’ 

29 

66 

32 

32 

37 

38 

32 

70 

37 

35 

42 

40 

%/2 

449 

584 

440 

450 

459 

458 

513 

638 

500 

515 

517 

525 

^Hnn 

114 

147 

105 

109 

107 

103 

136 

150 

125 

130 

119 

127 

^^7/2'’ 

306 

402 

319 

300 

294 

331 

344 

431 

359 

336 

373 

329 

®The  ^f5/2,  ^Hg/2,  and  ^f7/2  are  intermixed,  as  are  *05/2  and  ^Gy/2- 
‘’The  ^$3/2  multiplet  is  very  strongly  mixed  with  the  V7/2  multiplet. 

The  ^Gj/2  multiplet  is  strongly  mixed  with  the  ^G5/2  multiplet,  and  levels  belonging  to  ^G7/2  were 
arbitrarily  selected. 

Table  8.  Total  crystal-field  splittings,  AEj  (cm~^),  of  each  /  multiplet  of  Nd^"*"  in  SPAP  in  site 
with  =  -1.8  and  -2.“ 


qo  =  -1-8  qo  =  -2 


[(L,S)n 

A 

B 

Cl 

C2 

D1 

D2 

A 

B 

Cl 

C2 

D1 

D2 

%/2 

464 

586 

461 

451 

462 

449 

527 

619 

525 

512 

526 

505 

%l/2 

362 

493 

360 

340 

360 

321 

422 

511 

420 

400 

420 

381 

%3/2 

427 

580 

424 

405 

424 

385 

497 

605 

495 

475 

493 

455 

%5/2 

623 

807 

623 

600 

624 

580 

727 

1859 

727 

705 

727 

683 

%/2 

173 

309 

178 

178 

192 

195 

207 

347 

215 

211 

231 

226 

"53/2^= 

36 

49 

36 

37 

26 

27 

28 

53 

29 

30 

30 

33 

‘^Bg/2 

314 

457 

318 

317 

332 

336 

364 

495 

367 

362 

380 

378 

^Hnn 

81 

125 

78 

77 

77 

79 

96 

128 

93 

90 

92 

90 

2G7/2^ 

235 

369 

234 

246 

234 

258 

229 

393 

228 

265 

228 

250 

“The  ‘^F5/2,  ^H9/2,  and  ^¥7/2  are  intermixed,  as  are  '*G5/2  and  ^Gy/2- 
’’The  '^53/2  multiplet  is  very  strongly  mixed  with  the  *¥7/2  multiplet. 

'The  ^G7/2  multiplet  is  strongly  mixed  with  the  ^G5/2  multiplet,  and  levels  belonging  to  ^G7/2  were 
arbitrarily  selected. 


As  a  final  computation,  the  crystal-field  components,  were  calculated 
for  the  Ml  site  with  C3  symmetry  in  CPAP  and  SPAP,  and  the  crystal-field 
parameters  were  calculated  for  Nd^"’’.  In  these  calculations,  charge  com¬ 
pensation  was  ignored.  The  results  for  both  CPAP  and  SPAP  are  given  in 
table  9  along  with  the  invariants  I„(B)  for  each  case.  This  calculation  corre¬ 
sponds  to  type  A  in  table  4  for  the  M2  site.  On  comparing  I„(B)  from  table 
5  for  type  A  with  the  1„(B)  of  table  9  for  CPAP,  we  find  that  12(B)  and  14(B) 
for  Cj  symmetry  are  larger  than  12(B)  and  14(B)  for  C3  S3m[imetry.  However, 
the  opposite  is  true  for  the  corresponding  values  of  I^(B),  where  I^(B)  for 
the  Ml  site  in  CPAP  is  larger  than  /5(B)  for  the  M2  site  in  CPAP.  This  same 
result  is  true  for  the  I„(B)  in  SPAP.  The  energy  widths  of  the  multiplets  are 
given  in  table  10,  and  a  comparison  of  the  results  for  case  A  in  table  7 
shows  mixed  results;  that  is,  some  of  the  widths  of  the  multiplets  for  the 
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Table  9.  Nd^'*’  crystal- 
field  parameters, 

(cm^^),  and  rotational  - 

CPAP 

SPAP 

B„m,  ln(B) 

qo  =  -1-8 

n 

0 

00 

r-H 

II 

0 

qo  =  -'^ 

invariants,  I^{B)  (cm~^). 

B20 

1247 

1556 

936 

1173 

for  C3  (Ml)  site  in 

B40 

-2252 

-2359 

-1585 

-1646 

CPAP  and  SPAP  with 

B43 

1016 

1114 

700 

759 

qp  =  -1.8  and  -2. 

Beo 

-1176 

-1312 

-814 

-909 

890 

935 

629 

673 

ImBgs 

-260 

-432 

65.7 

194 

ReB66 

-117 

-3.61 

-130 

-35.9 

ImBee 

356 

417 

-249 

-311 

h(B) 

1247 

1556 

936 

1173 

14(B) 

2671 

2837 

1869 

1965 

h(B) 

1544 

2047 

1064 

1415 

1  able  10.  Total 
crystal-field  splitting, 
AEj  (cm“^),  of  each  / 

CPAP 

SPAP 

l(l,S)n 

00 

1 

n 

0 

11 

0 

qo  =  -1.8 

qo--'2- 

multiplet  of  Nd^'*'  in 

%/2 

667 

735 

471 

529 

CPAP  and  SPAP  in 

373 

426 

260 

298 

C3  (Ml)  site  with  qp  = 

%3/2 

439 

4892 

308 

347 

-1.8  and  -2. 

*h5/2 

861 

944 

595 

652 

^£3/2 

144 

167 

117 

139 

^S3/2 

26 

30 

17 

21 

%/2 

367 

422 

268 

308 

^Hn/2 

142 

165 

98 

114 

2G7/2 

324 

362 

233 

255 

Ml  sites  in  table  10  are  larger  than  the  corresponding  multiplets  for  the  M2 
sites  given  in  table  7.  The  results  for  SPAP  given  in  table  10  and  for  case  A 
of  table  8  are  similar  to  the  result  for  CPAP. 


4.  Discussion  and  Conclusion 

A  number  of  the  calculations  presented  here  were  made  in  an  attempt  to 
find  some  general  trend  and  eventually  a  simple  criterion  that  could  be 
used  to  sort  out  the  various  sites  in  the  experimental  data.  The  rotational 
invariants  are  one  such  quantity  that  I  chose  in  the  attempt  to  correlate  the 
^n(^)  with  the  AEj  of  various  multiplets.  In  particular,  if  we  select  the  ji 
splittings  and  compare  them  with  the  /2(B)  for  the  different  charge  com¬ 
pensations,  the  results  for  the  splittings  A£  (^£3/2)  are  B  >  D1  >  D2  >  C2  > 
Cl  >  A;  for  12(B),  on  the  other  hand,  B  >  A  >  C2  >  D1  >  Cl  >  D2.  If  the  crys¬ 
tal  field  were  weaker,  one  would  expect  a  one-to-one  correlation  between 
12(B)  and  AE  (‘^£3/2),  since  for  weak  crystal  fields,  A£  (“^£3/2)  must  be  pro¬ 
portional  to  12(B)  when  /-mixing  by  the  crystal  field  is  ignored.  However, 
that  /-mixing  cannot  be  ignored  is  obvious  from  the  above  results.  An 
examination  of  the  composition  of  the  ^£3/2  eigenfunctions  reveals  that  the 
^£3/2  lowest  energy  level  wavefunction  for  case  B  in  CPAP  (ijq  =  -2)  is 
approximately 
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90%  4f3/2  +  6%  4F5/2  + 1%  4G5/2 ; 
for  the  same  case,  the  highest  level  for  the  eigenfunction  is 

91%  4F3/2  +  3%  %/2  +  3%  %i2 . 

For  case  A,  the  wavefunctions  for  the  corresponding  levels  are 

95%  4F3/2  +  3%  4F5/2  +  0.50% 
for  the  lowest  energy,  and 

92%  4F3/2  +  5%  4F5/2  +  1%  %i2 

for  the  highest.  This  admixture  of  the  wavefunctions  of  different  /  values, 
caused  by  the  strength  of  the  crystal  field,  makes  a  simple  first  analysis  of 
the  experimental  observation  of  the  splitting  of  the  impossible. 
Obviously,  for  the  higher  energy  multiplets,  which  are  not  so  well  isolated 
from  one  another,  the  mixing  is  confusingly  severe:  For  case  A  in  CPAP 
{qo  =  -2),  the  wavefunction  for  the  lower  “^53/2  is 

71%  453/2  +  24%  %/2  +  2%  4F5/2 
and  for  the  highest  is 

83%  453/2  +  15%  4F7/2  +  1%  4F9/2 . 

Thus,  the  operator  equivalent  method  or  any  other  method  that  ignores  J- 
mixing  by  the  crystal  field  would  fall  far  short  of  the  mark  of  a  good  repre¬ 
sentation  of  the  experimental  situation. 

The  relationship  of  the  ^Ij  AEj  to  the  I„{B)  is  quite  complicated  even  in 
theory,  and  one  might  analyze  a  direct  relation  to  the  themselves  by 
using  some  of  the  results  given  by  Karayianis  [20]  where  /-mixing  by  the 
crystal  field  is  included. 

Finally,  a  calculation  of  the  crystal-field  parameters  for  Nd^'*'  in  the  Ml  site 
(table  1)  with  C3  symmetry  for  qQ  =  -1.8  and  -2  was  performed  for  the  case 
of  no  charge  compensation;  these  results  are  given  in  table  9  for  CPAP  and 
SPAP.  This  calculation  corresponds  to  the  A  type  of  calculation  of  table  4. 
The  f„{B)  values  in  comparison  with  the  corresponding  values  of  f„(B)  for 
the  Ml  in  table  9  show  that  l2iB)  and  /4(B)  are  larger  for  the  M2  site  than  for 
the  Ml  site,  but  the  reverse  is  true  for  the  corresponding  values  of  /^(B). 
However,  the  AEj  for  the  different  multiplets  of  the  4/  term  show  erratic 
behavior;  that  is,  for  the  qQ  =  -1.8  in  CPAP,  the  AE9 /2  for  M2  are  smaller 
than  those  for  Ml;  for  AE14/2,  M2  >  Ml;  for  AE43/2,  M2  >  Ml;  for  AE45/2, 
M2  <  Ml.  Similar  results  are  obtained  for  SPAP.  Thus  it  is  apparently  not 
possible  to  separate  the  experimental  data  by  sites  by  identifying  the  width 
of  the  splitting  of  even  well-defined  multiplets  unless  a  careful  analy¬ 
sis  of  these  splittings  using  some  of  the  results  of  Karayianis  [20]  revealed 
a  somewhat  simple  relationship  of  these  splittings  of  the  multiplets  with 
the  crystal-field  parameters. 
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